Resonant X-ray Scattering Measurements of a Spatial Modulation of the Cu 3d and 
O 2p Energies in Stripe-Ordered Cuprate Superconductors 
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A prevailing description of the stripe phase in underdoped cuprate superconductors is that the 
charge carriers (holes) phase segregate on a microscopic scale into hole rich and hole poor regions. 
We report resonant elastic x-ray scattering measurements of stripe-ordered Lai.475Ndo.4Sro.i25Cu04 
at the Cu L and O K absorption edges that identify an additional feature of stripe order. Analysis 
of the energy dependence of the scattering intensity reveals that the dominant signature of the 
stripe order is a spatial modulation in the energies of Cu 3d and O 2p states rather than the 
large modulation of the charge density (valence) envisioned in the common stripe paradigm. These 
energy shifts are interpreted as a spatial modulation of the electronic structure and may point to a 
valence-bond-solid interpretation of the stripe phase. 

PACS numbers: 74.72. Gh, 61.05.cp, 71.45.Lr, 78. 70. Dm 



Static stripe order in cuprates was first theoreti- 
cally predicted by mean-field Hubbard model calcula- 
tions [IHl] and subsequently observed in lanthanum- 
based cuprates by neutron and x-ray diffraction [SHU]. 
Although still a matter of debate, more recent work has 
indicated that stripe-like density wave order is generic to 
the cuprates [12H17' and plays a significant role in com- 
peting with or possibly causing superconductivity 18 . 

Microscopically, stripes in the cuprates have been 
widely described as rivers of charge — hole-rich antiphase 
domain walls that separate undoped antiferromagnetic 
regions. However, alternate models with different under- 
lying physics, such as the valence bond solid (VBS), have 
also been proposed to explain stripe order [I9 H 2T ] . VBS 
models involve singlet formation between neighbouring 
spins and, in contrast to other models of stripe order, 
may occur with a small modulation of the charge density 

Distinguishing which of these models is most relevant 
to stripe order in the cuprates is challenging since the 
models share many symmetries and experimental signa- 
tures. In particular, direct evidence for charge-density 
modulations, which may distinguish various models, has 
been elusive. Neutron and conventional x-ray scattering 
are only sensitive to lattice displacements. It is therefore 
only inferred indirectly that these lattice displacements 
are induced by modulations in charge density (valence). 
Resonant soft x-ray scattering (RSXS) offers a means 
to couple more directly to modulations in the electronic 
structure, including charge density modulations. By per- 
forming an x-ray diffraction measurement on resonance 
(at an x-ray absorption edge), the atomic scattering form 
factor /(w) is enhanced and made sensitive to the va- 
lence, orbital orientation and spin state of specific ele- 
ments. A key feature of RSXS is that the energy depen- 



dence of the scattering intensity through an absorption 
edge differs for lattice distortions, charge-density modu- 
lations or other forms of electronic ordering, providing a 
means to distinguish these different types of order. 

In the cuprates, RSXS of the [2e, 0, L] charge 
density wave (CDW) superlattice peak has been 
measured in stripe-ordered La2_j;Baa;Cu04 (LBCO) 
[lO], La2-:,-yEuySr^Gu04 (LESCO)[ni El and 
Lai.475Ndo.4Sro.i25Cu04 (LNSCO) ^ at the O K 
(Is — >■ 2p) and Cu L (2p —>■ 3d) absorption edges, which 
provide sensitivity to the O 2p and Cu 3d orbitals 
that are central to the physics of the cuprates. These 
measurements have been interpreted as direct evidence 
for a large valence modulation on the O sites jlO) . 
Moreover, it is argued that a modulation of the valence 
occurs primarily on the O sites and not on the Cu sites, 
which are instead subject to lattice distortions induced 
by the valence modulation on the O sites [101 fTT] . 
However, efforts to model the energy dependence of the 
scattering intensity based on this picture are not truly 
reconciled with experiment, leaving this interpretation 
open to question [TT] . 

In this Letter, we present O K and Cu L edge RSXS 
measurements of LNSCO. The energy dependence of the 
scattering intensity is modelled using x-ray absorption 
measurements to determine /(w) at different sites in the 
lattice, a procedure that has proven effective in describ- 
ing the scattering intensity of valence modulations in the 
chain layer of ortho-IIYBazCugOe+a (YBCO) gl|- Con- 
trary to previous analysis of LESCO [11] and LBCO [lO] . 
we show that the resonant scattering intensity is best de- 
scribed by small energy shifts in the O 2p and Cu 3d 
states at different Cu and O sites rather than a valence 
modulation of O and a lattice displacement of Cu. 

The measured intensity for scans through Q — [H, 0, 
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FIG. 1. (color online). Resonant scattering of 1/8 doped LNSCO at the Cu L and O K edges, (a) Schematic representations 
of bond-centred stripe ordering for the three models proposed to describe the resonant scattering energy dependence, (b) H 
scan through the CDW superlattice peak at [H, 0, 1.5] and at the peak of the Cu L3 absorption edge |24j. (c) — (h) Scattering 
intensity as a function of photon energy through the Cu L3,2 and O K absorption edges. The measured intensity (red line) 
is compared to the scattering intensity of valence modulation (blue line), lattice displacement (green line) and energy shift 
(black line) models of the stripe-ordering. The best agreement with experiment is for the energy shift model, (i) and (j) X-ray 
absorption fi{uj) at the Cu Z/3,2 (i) and O K (j) absorption edges measured using total electron yield. 



1.5] vifith the photon energy at the peak of the Cu L3 
absorption edge (931.3 eV) are presented in Fig. 1(b) 
(see Supplemental Material [5S] for methods). Below 
the stripe-ordering transition temperature of ^70 K, a 
clear superlattice peak is observed at Q = [2£,0, L] = 
[0.236,0,1.5]. Above the stripe-ordering transition tem- 
perature, a large, smoothly varying fluorescence back- 
ground is observed. To determine the photon energy de- 
pendence of the scattering intensity, H scans through 
the superlattice peak were performed at 22 K as a func- 
tion of photon energy at L = 1.5 for Cu and L = 0.2 
for O. The fluorescence background is fit and subtracted 
from each scan. The integrated intensity of the super- 
lattice peak is then determined as a function of photon 
energy for the O K and Cu L edges [Figs. 1(c) and 1(d)]. 
The resulting spectra are qualitatively similar to previous 
measurements on LBCO y^ and LESCO [11]. Impor- 
tantly, our measurements extend the previous Cu L edge 
measurements to include the L2 edge, which proves valu- 
able in distinguishing models for the stripe phase. An 
important feature of our measurement is that all scatter- 
ing measurements are performed with the incident x-ray 
polarization along the bo axis of the sample. As a result, 
the scattering intensity will be sensitive to only the O 
Py and not the O Px orbitals. Assuming doped holes go 
only into CT-bonded orbitals of O, this measurement ge- 



ometry is only sensitive to half of the oxygen atoms; the 
site-centred and not the bond-centred oxygen. This fact 
simplifies the expression for the structure factor. 

The measured energy dependence of the scattering 
intensity is compared to three model calculations [see 
Fig. [l[a)]: 1. valence modulation, a spatial modulation 
in the valence of the Cu and O; 2. lattice displacement, 
a small displacement of the Cu and O atoms from their 
equilibrium positions outside the stripe-ordered phase; 
and 3. energy shift, a spatial modulation in the energy 
of the Cu 3d and O 2p states. The first two models es- 
sentially follow previous analysis of RSXS in LBCO and 
LESCO [lOllII]. 

The three models differ in the structure factor (de- 
scribed in the Supplemental Material [5S]) and the energy 
dependence of the atomic scattering form factor /(w). 
These two factors give rise to a different energy depen- 
dence to the scattering intensity, /^(w). For the valence 
modulation model, Is{uj) oc \f{uj,p2) — /(w,pi)|^//i(a;), 
where pi and p2 are the local hole concentrations (va- 
lence) at different sites [see Fig. [IJa)] and ^(w) is the ab- 
sorption coefficient. For the lattice displacement model, 
/(w) is the same at each site for a given element and 
/s(a;) ex |/(cj)|^//i(a;). Finally, for the energy shift model, 
Isioj) cx \f{hu; + AE) - f{huj - AE)\'^ / fi{uj) , similar to 
the valence modulation model but with an energy shift 



zLAE at different sites instead of a modulation in va- 
lence. In all three models, the site specific f{uj,pi^2, ^E) 
are determined from x-ray absorption measurements. 

Valence modulation model. — For the valence modula- 
tion model, x-ray absorption spectroscopy (XAS) on sam- 
ples with different doping levels from Ref. |26j are used 
to determine f{uj,pi^2)- This procedure found very good 
agreement between experiment and calculations for a 
modulation of the Cu valence in the chain layer of oxygen- 
ordered YBCO [21]. In lanthanum based cuprates, the 
key features of the O K edge XAS are two preedge peaks 
at 528.6 eV and 530.5 eV that are due to hybridization 
between Cu 3d and O 2p states and have been assigned 
to the mobile doped holes and the upper Hubbard band 
respectively [Fig. [T]^j)] [26H28] . The intensities of these 
two peaks evolve strongly with doping, whereas the spec- 
tra at higher energy are doping independent and domi- 
nated by O 2p states hybridized with rare earth 5d and 4/ 
states [25J . As argued in Ref. pT| , the scattering intensity 
for a valence modulation of arbitrary magnitude can be 
modeled using XAS measured at two different dopings. 
Here fj{u!) [Fig. 2(b)] and the scattering intensity ex- 
pected for a valence modulation [Fig. 1(d)] is calculated 
from XAS in La2_:rSr^Cu04 (LSCO) at a; = 0.07 and 
0.15 from Ref. [26 , corresponding to a hole modulation 
of Sp — pi — p2 — 0.08. Although this calculation suc- 
cessfully produces two peaks at approximately the cor- 
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FIG. 2. (color online). The atomic scattering form factors 
f{io) as a function of photon energy through the Cu L and O 
K absorption edges for (a), (b) La2-2!Sri:Cu04 {x — 0.07 and 
0.15) and (c), (d) Lai.475Ndo.4Sro.i25Cu04. f{uj) for LSCO 
is determined using XAS measurements from Ref. [26] . 



rect energies, it strongly overestimates the intensity of 
the peak at 529.9 eV [TT]. A different choice of doping 
values to determine /(a;,pi_2) impacts the magnitude of 
the scattering, scaling it roughly as Sp^, and produces 
differences in the line shape (see Supplemental Material 
[25]). However, calculations using existing XAS data are 
all similarly inconsistent with the measured RSXS line 
shape. 

A similar analysis, again using XAS from Ref. [26] to 
determine f{uj,p) [Fig. 2(a)], can be applied to the Cu L 
edge. The XAS for the Cu L edge exhibits two primary 
peaks at 931.3 eV and 951.3 eV corresponding to the L3 
and L2 edges that are split by the spin-orbit coupling 
of the 2p core electrons. Focusing on the L3 edge, the 
XAS is comprised of a peak (931.3 eV) and a shoulder 
(932.3 eV) that are associated with dP (a single hole in 
the d^2_y2 orbital) and d^L (doped holes that are pri- 
marily on the oxygen ligands) ground states. Consistent 
with this assignment, the high energy shoulder evolves 
much more strongly with doping than the d^ peak [26| . 
It follows that the predicted scattering intensity for a 
valence modulation of the Cu is peaked at the shoulder 
and not the peak of the XAS [Fig. 1(c)]. As discussed in 
Refs. [ini [H] , this is in poor agreement with the energy 
dependence of the resonant scattering, which is peaked 
at the maximum of the XAS. 

Lattice displacement model. — Calculations based on 
the lattice displacement model [Figs. 1(e) and 1(f)] are 
also in poor agreement with experiment [here using XAS 
on our sample of LNSCO in Figs. l(i) and l(j) to deter- 
mine /(w) in Figs. 2(c) and 2(d)]. The lattice displace- 
ment model at the O K edge predicts large scattering in- 
tensity above and below the absorption edge that is not 
observed in experiment and, at the Cu L edge, scatter- 
ing intensity that is broader in energy and has a smaller 
ratio of the L3 to L2 peaks than the measurement. The 
calculated magnitude of the scattering intensity assumes 
a 0.004 A lattice displacement, as deduced from neutron 
scattering [5]. 

Energy shift model. — Like the lattice displacement 
model, XAS on our sample of LNSCO [Figs. l(i) and 
l(j)] is used to determine f{uj) [Figs. 2(c) and 2(d)] for 
the energy shift model. For the two sites (1 and 2), /(w) 
is shifted in energy by AE = ±0.1 eV for both the O K 
and Cu L edges. In contrast to the lattice displacement 
and valence modulation models, the energy shift model 
is in very good agreement with experiment. At the Cu L 
edge, it captures the correct intensity ratio of the Cu L3 
and L2 peaks, the correct width in energy of the scatter- 
ing and the correct energy position of the maximum of 
the scattering intensity. Similarly, at the O K edge, the 
energy shift model reproduces the energy dependence of 
the preedge peak. It does not agree with the spectra at 
higher energy, predicting a large peak at 531.5 eV that 
is not observed. However, this apparent discrepancy can 
be reconciled if we interpret this as evidence that only 



the low energy states involving hybridization between 
the O 2p and Cu id states (and not the rare earth 5d 
and 4/ states) are subject to these energy shifts. We 
also note that the choice of A_E impacts the magnitude, 
which scales as Ai5^, but not the energy dependence of 
the calculated scattering intensity, provided AE is less 
than the energy width of the XAS ('-^0.2 eV). As such, 
Ai? is neither determined in our analysis, nor should it 
be viewed as a fitting parameter. 

Our interpretation of the energy shifts is that they are 
induced by subtle spatial modulations of the local elec- 
tronic structure. The energy levels of the unoccupied 
Cu 3d and O 2p states can be described by parameters 
such as the Cu onsite Coulomb repulsion [Udd), the crys- 
tal field splitting parameters {10Dq,Ds,Dt), the charge 
transfer energy (Apd) and the Cu 3d — O 2p hopping {tpd) 
P71 1501 l5Tj . Small changes to these parameters can lead 
predominately to shifts in the Cu 3d and O 2p energy 
levels that would manifest themselves as energy shifts in 
/(w) j32j, unlike the large changes in the spectral weight 
distribution that are observed with doping |26| . Since the 
XAS and RSXS are sensitive to the Cu 2p —>■ 3d and O 
Is— ^2p transitions, modulations of the Cu 2p and O Is 
binding energies may also contribute to the energy shifts. 

The agreement between the measured scattering inten- 
sity and this simple phenomenological model indicates 
that energy shifts are responsible for the dominant con- 
tribution to the resonant scattering intensity. In compar- 
ison, contributions arising from lattice displacements and 
valence modulations appear to be much less significant. 
This is reasonable for the lattice displacement model, 
given that the magnitude of the calculated scattering in- 
tensity is ^4 orders of magnitude weaker than the energy 
shift model at both the O K and Cu L edges [Since AE 
is unknown, for this analysis we compare to an arbitrar- 
ily chosen value of AE — 100 meV, which serves as a 
reasonable upper limit value. Comparing to other values 
of A^ involves scahng by (A£'/(100 meV))^] 

In contrast, the valence modulation model (assuming 
6p = 0.08) is ~4 times larger than the energy shift model 
at the O K edge. At the Cu L edge, the valence modu- 
lation model is ~35 times weaker than the energy shift 
model at their respective peak values but is comparable 
in intensity at 932.3 eV, the peak energy of the valence 
model calculation. As such, unlike lattice displacements, 
we do not expect the energy shift contribution to dom- 
inate the scattering intensity for valence modulations of 
order Sp = 0.08 or larger. This argues against a large va- 
lence modulation, such as those reported in Ref . [10] , but 
does not rule out smaller valence modulations. However, 
even if negligibly small, one can infer that valence modu- 
lations must be nonzero, as they must occur for a spatial 
modulation of the Cu 3d and O 2p energies (i.e., the en- 
ergy shifts provide indirect evidence for valence modula- 
tions) . Placing a precise upper limit on the magnitude of 
the valence modulation is beyond the scope of the present 



work, requiring more sophisticated modeling. However, 
we note that the energy shift model has the same unoccu- 
pied spectral weight, and hence the same valence, for all 
sites in the stripe phase. As such, the resonant scattering 
line shapes, which are well described by the energy shift 
model alone, are consistent with a stripe phase that has 
a minimal but nonzero valence modulation. 

The origin of these modulating energy levels and how 
they relate to the microscopic mechanism of stripe order- 
ing is an open question. The energy shifts may simply 
be induced by small charge-density modulations or lat- 
tice displacements, yet still be the dominant signature in 
resonant scattering. Alternatively, they may be a more 
direct signature of the interactions underlying stripe or- 
der. For instance, these energy modulations may point 
to a VBS description of the stripe phase [T^21j . In the 
VBS picture, stripe order is driven by exchange interac- 
tions, which also induce lattice displacements and bond- 
centered charge order. However, the magnitude of bond- 
centered charge density modulations can be small, being 
screened by long-range Coulomb repulsions. This may 
provide an explanation for the lack of clear evidence for 
a valence modulation from resonant scattering. In con- 
trast, the energy shifts that we have identified in RSXS 
may arise naturally out of a VBS, which involves a mod- 
ulation of the bonding in the lattice. 

Finally, our identification of energy shifts is likely ap- 
plicable not only to CDW order in other cuprates (en- 
ergy shifts were recently shown to also describe den- 
sity wave order in YBCO [53]), but also to other transi- 
tion oxides. For instance, recent first principles calcula- 
tions have shown that several "charge" ordered transition 
metal oxides exhibit a site dependence to the energies of 
the electronic states but no site dependence to the total d 
orbital occupation [34 , similar to our phenomenological 
energy shift model. 
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